Monitoring Quantum Oscillations with very small Disturbance 
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We present a new scheme to detect and visualize oscillations of a single quantum system in real 
time. The scheme is based upon a sequence of very weak generalized measurements, distinguished 
by their low disturbance and low information gain. Accumulating the information from the single 
measurements by means of an appropriate Bayesian Estimator, the actual oscillations can be moni- 
tored nevertheless with high accuracy and low disturbance. For this purpose only the minimum and 
the maximum expected oscillation frequency need to be known. The accumulation of information 
is based on a general derivation of the optimal estimator of the expectation value of a hermitian 
observable for a sequence of measurements. At any time it takes into account all the preceding 
measurement results. 



I. INTRODUCTION 

Consider a two-level quantum system where the prob- 
ability to find the system in a projection measurement 
on a specific level oscillates due to a periodically time- 
dependent external potential. We propose a measure- 
ment scheme which allows to monitor these oscillations 
by means of a sequence of consecutive measurements car- 
ried out on a single two-level system. As in the case of 
the detection of gravitational waves we assume that we 
arc dealing with a one shot experiment, i.e. the mea- 
surements have to be carried out on a single quantum 
system and the experiment cannot be repeated in order 
to acquire more measurement data. The proposed mea- 
surement scheme yields a real time record of the actual 
oscillations, which are sometimes called Rabi oscillations. 
This is done by appropriately estimating after each mea- 
surement in the sequence the actual value of the oscillat- 
ing probability. 

In contrast to measurements in classical physics, quan- 
tum measurements have the following prominent fea- 
ture: The more information they provide, the more they 
change the state of the measured system. The most pre- 
cise measurements are von Neumann projection measure- 
ments. They project a quantum system in an eigenstate 
of the measured observable, which represents in general 
a drastic disturbance of the system's state. 

There is a broader class of measurements called gener- 
alized measurements, which can be realized by coupling 
the system in question to an ancilla system and carrying 
out a projection measurement on the ancilla. Depend- 
ing on the kind of coupling and its strength these indi- 
rect measurements can exert an influence on the system 
which ranges from very weak to very strong. They can 
be elegantly described in the POVM formalism 0]. In 
order to keep the disturbance caused by a sequence of 
measurements low, we employ generalized measurements 
with very weak influence. On the other hand, the weaker 
the influence of these measurements is, the less infor- 
mation about the measured system they convey. This 
disadvantage can be compensated by accumulating data 



from all measurements. A corresponding data processing 
scheme is proposed below. 

The strength of the influence of a sequence of measure- 
ments can be determined by considering the case where, 
apart from the measurements, no other dynamical influ- 
ence is present. In the example of the oscillating two-level 
system this corresponds to screening or turning off the 
time-dependent potential. For an appropriate sequence 
of measurements a measure for the strength of the in- 
fluence is then given by the "decoherence time". That 
is the period after which the coherences of the systems' 
state have decayed to 1/e of their initial value (for qubits 
cp. 0). The undisturbed dynamics of the system, on the 
other hand, can be characterized by the time scale Tr, 
which is the period of the oscillations, if no measurements 
are carried out. 

Comparing the decoherence time Td to the period Tr, 
roughly three modes or regimes of measurement can be 
distinguished Q : (i) for Td ^> Tr the system evolves ap- 
proximately according to its undisturbed dynamics, i.e. 
the disturbing influence of the measurements is compara- 
tively small, (ii) Td ~ Tr, both dynamical influences are 
equally strong, (iii) Td -C Tr, decoherence induced by 
the measurements dominates the dynamics of the system. 
Looking in mode (iii) at the systems dynamics in the se- 
lective regime, i.e. given certain measurement outcomes, 
one finds quantum jumps or in the limit of a continuous 
projection measurement the Quantum Zcno effect, where 
the system freezes in an eigenstate of the measured ob- 
servable. 

In it was shown that a detection of the oscillations 
of a two-level system under the influence of an external 
field with reasonable accuracy and disturbance is possi- 
ble employing mode (ii). A physical realization of this 
measurement scheme was proposed by probing a photon 
oscillating between two cavities by a sequence of Rydberg 
atoms 0, based on an experiment of Haroche et al. ||. 

In contrast to the latter investigations we want to show 
that the results can be improved by working in mode 
(i). The advantage of mode (i) obviously is the weak 
influence and thus the low disturbance inflicted by the 
measurements. On the other hand mode (i) represents a 
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challenge because there the measurement results are only 
poorly correlated to the actual state evolution. In order 
to overcome this difficulty we study optimal estimates 
in sections ITfl and ITTT1 The presented approach is rather 
general. We optimally estimate the expectation value of 
an arbitrary hermitian observable of a quantum system 
with finite dimensional Hilbert space, given the result of 
a generalized measurement (which can also consist of a 
sequence of consecutive measurements). The result can 
be applied to our special case to estimate the probability 
to find a two-level system in a projection measurement on 
a specific level. This is done in section llVl after giving a 
brief description of this scheme. The results of numerical 
simulations of our measurement scheme are discussed. 
An appendix contains a recipe for these simulations with 
useful formulae to abbreviate the computations. 



II. ESTIMATOR FOR MEAN VALUES OF 
OBSERVABLES 

There probably exists considerably more literature on 
state-estimation than on the estimation of mean values 
of physical quantities such as energy, position or spin 
of quantum systems. Nevertheless there might be ques- 
tions which do not require the maximal knowledge of the 
statistics of all measurements that can be carried out on 
a quantum system — as it is represented by the state of 
the system — but rather the knowledge of a single phys- 
ical property such as the mean position of a quantum 
particle. For example, when attempting to detect gravi- 
tational waves, only the mean spacial distances between 
test masses have to be estimated at different times @ . Of 
course the calculation of the mean value of an observable 
as well as state determination is not an issue if a large en- 
semble of identically prepared systems is available to be 
measured. But for experiments with restricted resources 
and especially in one-shot experiments estimation proce- 
dures become essential. 

Let us consider the following task. Given a quan- 
tum system with d-dimensional Hilbert space Ti. After 
a POVM measurement with result m the state of the 
system reads 



Mr, 



(i) 



where M rn is the Kraus operator corresponding to the 
measurement result m, and 



p(m|V) = (MMlM m \^) 



(2) 



represents the probability to obtain m, provided the state 
before the measurement was Let us assume that we 
know the result m and the respective Kraus operator M m 
but we do not know the initial state . What is the best 
way to estimate the parameter 



which represents the expectation value of the observable 
A with respect to the state l^m)? 

It turns out to be rewarding to base the parameter es- 
timation on the least squared error criterion: The value 
g m is an optimal estimate of the parameter 9 m if it min- 
imizes the expected square of the error 



E ((9 m — g n 



J (fim - .9m) 2 p(^|V')p('0)rf'0 
J p(m\t/})p(ijj)dijj 



(4) 



For the sake of simplicity we assume no prior knowledge 
about the initial state \ip), i.e. p(ip) = 1 and dip is a 
normed measure over the set of all pure states which is 
invariant under the action of the rotation group SU (d) . 

Taking into account the linearity of the expectation 
value, it is easy to see that in this case the optimal esti- 
mator g m is equal to the expected value of 9 m : 



E((6 



g m f) = E ((9 m ) 2 ) - 2gE(6 m ) + g 2 m 

,2 . ,r ,„ x ( 5 ) 



(E(9 n 



Var(# m ) , 



where 



Y^{e m ) = E{e 2 m ) - E{e v 



(6) 



represents the variance of 9 m . The right-hand side of (JSJ 
assumes a minimum for g m = E(9 m ). Such a value g m is 
also called the Bayesian estimate (cp. Q). 
Evaluating g m we obtain: 

_ J 9 m p(m\ip)dip 
J p(m\tp)di() 

J (ip\MrliM m \ip}dtp 
= tr[MlAM m } 
tr[AGt»Af m ] 



(7) 



{ip m \A\ijj ri 



(3) 



This quantity represents the best estimate of the expec- 
tation value of observable A after one single generalized 
measurement with result m, if the state before the mea- 
surement is completely unknown. It is the best estimate 
in the sense that it leads to the least expected squared 
error. 

Formula (JJJ can in particular be applied to estimate, 
after a generalized measurement, the probability to find 
a system with two levels and 1 on level 1 . In this case A 
should be chosen to be the projector on level 1 since the 
expectation value of this projector is equal to the desired 
probability. 



III. ESTIMATOR FOR SEQUENTIAL 
MEASUREMENTS 



In this section we derive an optimal estimator for a 
sequence of measurements. For the sake of broad ap- 
plicability we consider the general case of a sequence of 
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generalized measurements carried out on a single d-level 
system with unknown Hamiltonian. 

The single measurements with Kraus operators N mn 
are carried out consecutively on a single quantum sys- 
tem at times t = ht, where n is an integer. Here the 
number m n represents the result of the n-th measure- 
ment. The single measurements are of duration St and 
during this time the motion due to the system's Hamil- 
tonian H can be neglected (impulsive measurement ap- 
proximation). Between two consecutive measurements, 
the system evolves according to the unitary operator 



U = cxp —Ht . 
h 



(8) 



For later convenience we express the unitary evolution 
by means of a unit vector k g M. d _1 with components kj 
and an angle < <f> < 2tt: 



U 



exp —Ht = cxp ik. ■ eq 
n 



(9) 



where k • e = y\ kjej, and the ej form a complete set of 
generators of SU(d). The state \i/j m ) of the system after 
n measurements with results (mi, . . . ,m n ) = m is then 
given by equation Q with Kraus operator 



M m = N mn UN mn _ 1 U, 



,N mi U 



(10) 



instead of M m . 

As above we assume that we know the results of the 
n measurements and, consequently, the corresponding 
Kraus operators N mi , . . . , N mn , but we know neither the 
unitary evolution U between the measurements nor the 
initial state of the system. Our ignorance about U 
has to be incorporated into the optimal estimate of the 
observable A after the n-th measurement. Instead of only 
averaging over all possible initial states as in Q we also 
have to average over all possible unit vectors k and possi- 
ble angles <j> weighted by the corresponding probabilities 
f>(k) and p(4>): 

= I dmP(m\iP, k, 4>)p{K)p{4>)dH d2 - 1 ^ 
9m Jp(m|V',k,0)p(k)p((/))dVd d2 - 1 k# ' 

We assume that the direction of k and the angle cf> are 
equally distributed, i.e. p(k)=const and p((f>)= const. 
The optimal estimator g m is then given by 



/ 9 m p(m\tp,k, <t>)dij)d d _1 k# 
J p(m\ip, k, <j))d'ipd d2 ~ 1 'kd(l) 

[ trjM^AM^d^^kdcj) 
J tx[MLM Ia ]d< p - x 'kd^ 



(12) 



For unknown unitary evolution, g m represents the esti- 
mate of the expectation value of observable A after n 
measurements with results m = (mi, . . . ,m n ). It mini- 
mizes the expected squared error. Calculating g m after 
each measurement in a sequence of measurements yields 
an optimally updated estimate of the current mean value 
of observable A. 



IV. APPLICATION: TRACKING AN 
OSCILLATING QUBIT 

We will now apply the estimator given in equation l|12(l 
to the sequential measurement of an oscillating qubit. 
Oscillating qubits arc realized for example by two-level 
systems such as coupled quantum dots, trapped atoms 
in an external field or photons oscillating between two 
microwave cavities Q. For the sake of concreteness, we 
consider a two-level atom under the influence of a reso- 
nant laser field. The Hamiltonian of such an atom can 
be approximated by 



H = £ o |0)(0| + £i|l)(l 



(13) 



(|1)(0| cxp{-iw<} + |0) (1| exp{iut}) , 

where uj = {E\ — Eo)/H, and the Rabi frequency flu 
represents the strength of the coupling between the atom 
and the electromagnetic field. The resulting motion of a 
state that is not subjected to measurement is represented 
by 



m)) =co(i)|0) + Cl (i)|l) 



(14) 



with | ci | 2 = |(1 + acos(S7fl.i + tp)) where the constants a 
and tp depend on the initial state of the qubit. 

In order to observe the actual behaviour of the system 
in real time, we have to measure and estimate the expec- 
tation value of |1)(1|, which is equal to |ci(t)| 2 . In the 
following we take the viewpoint that we already know the 
Hamiltonian (|13|l apart from the precise value of the cou- 
pling strength f2^j. This is quite a natural assumption, 
since the form of H has to be known in order to design 
measurements, which requires the knowledge of how to 
couple a meter to the system. 

In a first step we choose the Kraus operators. To ob- 
tain information about the observable |1)(1|, it seems 
natural to employ measurements the effects (N} n N m ) of 
which commute with 1 1)(1 1 (cp. ||). Because such effects 
are diagonal with respect to the basis states |0) and |1), 



NlNr, 



: P «|0)(0| 



(m) 
■Pi 



(15) 



the probability to obtain the result m depends directly 
on |ci(t)| 2 : 

pHv(f)) = (m\NiN rn \m), 

= p^ + Ap^laW 



(16) 



where Ap^ rn) := p\" L> —pi^ ■ For the sake of simplicity we 
consider measurements with two possible results, + and 

— , and Kraus operators N± = U±^J N±N± with trivial 
unitary part, i.e. U± = 1. The Kraus operators of each 
single measurement thus read 

N + := ^|0)(0| + yp+|l)(l|, (17) 

N- := xfp^\0)(0\ + Jp^\l)(l\ (18) 
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with positive numbers p- which satisfy Pj+Pj = 1 • 
A detailed analysis of optimal Kraus operators from the 
viewpoint of Bayesian estimates will be presented else- 
where. 

The change of state caused by consecutive measure- 
ments with Kraus operators N± as given in Eqs. (|17I18|) 
can be quantified by the decoherence time T<j. That is 
the average period after which the off-diagonal elements 
(i\p{t)\j) with i 7^ j G {0, 1} of the systems density oper- 
ator p(i) are decayed to 1/e of their original value, if the 
only dynamical influence is given by the measurements. 
T d is related to the decoherence rate 7 II by 



2 _ 8Tp(l-p) 
7 " (Ap) 2 



(19) 



withp := (j)q+Pi )/2 and Ap := p± — Pq ■ t is the time 
which passes between two consecutive measurements. If 
unitary dynamics generated by the Hamiltonian H <|13|) 
are present, they can create coherences representing a 
counter weight to the decoherence caused by the mea- 
surements. For 



T d » T R , 



(20) 



where Tr = 2^/17^ is the period of the Rabi oscillation, 
the influence of the single measurements on the state of 
the system becomes negligibly small as compared to the 
influence of the unitary dynamics. In this mode, which 
was called mode (i) in the introduction, we run our se- 
quential measurement. 

Apart from condition l|2U|l there is another require- 
ment for the sequence of measurements: a reasonably 
high number of measurements should take place on the 
time scale Tr of the unitary dynamics in order to resolve 
these dynamics. Hence, 



Tr » r . 



According to our experience based on numerical simula- 
tions Tr/t w O(10) is sufficient to obtain good results 
(sec below). Note that in order to meet conditions l|2(J|) 
and (|2T|) only a vague knowledge about the order of mag- 
nitude of the time scale Tr and accordingly of Qr is nec- 
essary. Knowing an upper bound and a lower bound 
Tft for Tr, both conditions can always be satisfied by first 
inserting the lower bound T^ into l|21|l and choosing t 
accordingly. Having fixed the value of r, the parameters 
Po and pi can be tuned such that 



8p(l-p)/(Ap) 2 »T>/r. 



(22) 



In other words: the sequential measurement has to con- 
sist of frequent measurements which are sufficiently weak. 

Since the influence of the measurements obeying H2U|) 
is very weak on the time scale Tr they convey only very 
little information about the system over the period Tr. 
This is where the accumulation of information by means 
of the Bayesian estimate (I12|) comes into play. The esti- 
mate after the n-th measurement makes the best possible 



use of the data collected in the previous n — 1 measure- 
ments according to the least square criterion. 

In our special case the Bayesian estimate given in (TT21) 
reduces to 



J tr[A/t,AM m ](ty 
/ tv[MLM m ]d<f> 



(23) 



with Kraus operators M m given by Eq. (|10|l . (|17(l and 
|T£)l. The unitary evolution U in M m (cp. Eq. flT))) ) is 

most easily represented in the interaction picture, where 
Hamiltonian l|13[l reads Hj = hQ,R(j x /2 with the Pauli 
spin operator a x := |1)(0| + |0)(1|: 



Z7 = exp(-i-|-) . 



(24) 



Here the angle <\> of rotation on the Bloch sphere is given 
by 



VLrt 



2kt 
~Tr~ 



(25) 



Because of condition l|21|l . which guarantees the tem- 
poral resolution of the Rabi oscillations by the sequence 
of measurements, the angle is in fact very small. We 
did not include this a priori information into the estimate 
g m used in our numerical simulations. Instead, we let <j> 
run from to 2tt in the integral in equation (|12|l . This 
corresponds to a completely unknown Rabi period Tr. 
An appropriate change of the range of integration might 
lead to an improvement of the estimate g m . 

Note that the representation of the Kraus operators 
N± in the interaction picture is the same as in the 
Schrodingcr picture used so far: 



N { ± n (t) = et^-^ATte-t^o (*-'<>) = N ± 



(26) 



(21) with H := E o \0}(0\ +E 1 \1)(1\. 



V. RESULTS AND CONCLUSION 

In Fig^the result from a numerical simulation of a se- 
quential measurement as specified above is plotted. The 
simulation started with the initial state \tp(to)) = |1). 
The solid curve represents the evolution of the estimate 
g and the dashed curve corresponds to the dynamics of 
|ci| 2 , taking into account the influence of the measure- 
ments. This represents what really happens. The dot- 
ted curve displays how |ci| 2 would have evolved without 
measurements . 

In the beginning of the sequential measurement (upper 
picture in Fig^l the curves of |ci| 2 with and without 
measurements are close; thus the disturbance due to the 
measurements is small. A small phase shift of the oscilla- 
tions however is recognizable. The values of the estimate 
g are not well correlated to the values of \c\ | 2 . After many 
measurements (middle picture in Fig.Q) the estimate g 
starts to approximate the disturbed |ci| 2 values, while 
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FIG. 1: Comparison between simulated evolution of a qubit's 
Rabi oscillations and processed measurement signal for p — 
0.5, Ap = 0.1 and r = T«/16. Dashed curve: |ci| 2 over 
time (in units of the Rabi period Tr) in the presence of weak 
measurements. Dotted curve: \c\ | 2 over time in the absence of 
measurements. The solid curve corresponds to the evolution 
of the estimate g based on the measurement results. 

the phase shift between the latter and the values of |ci| 2 
without measurements has increased. The amplitude of 
the oscillation of |ci| 2 is not changed by the measure- 
ments if — as seen here — it equals one in the absence of 
measurements. Eventually (lower picture of Fig^) es- 
timate g and the curve of actual values of |ci| 2 nearly 
coincide. 

The main results of the simulations are: i) With in- 
creasing time the estimate g reflects the actual oscilla- 
tions with growing fidelity. After approximately hundred 
Rabi cycles these oscillations are monitored by g with 
high accuracy, ii) In the presence of the weak measure- 
ments the sinusoidal shape and the period of the Rabi 
oscillation is almost the same as in the absence of mea- 
surements. The measurements however cause a phase 
shift of the oscillation. Therefore the estimate g also 
reflects shape and period of the undisturbed Rabi oscil- 
lation. 

This demonstrates that our approach allows to monitor 
the periodic evolution of an expectation value with high 
fidelity. The key to the monitoring are measurements 
with very low disturbance combined with an estimator 
which accumulates at any time the information gained in 
the sequence of previous measurements. 

VI. APPENDIX 

Simulations of the tracking procedure explained in sec- 
tion \lV\ were performed with a program ||| that is based 



on the following algorithm. 

1. Initialize the qubit's state vector \ip) and the num- 
ber of measurements, n max . Set n = 1. 

2. Evolve \ip) in time: \tp) -> e~ iHT / h \ip). 

3. Perform measurement: 

(a) Generate a (pseudo) random number m n 
whose value is either or 1, depending on the 
probability p mn = {^N^N^ \ tp). 

(b) Update \if>) -> {N mn I jK^M ■ 

4. Calculate the estimator according to formula J53J). 

5. If n < rimax; then continue at step |21 and increment 
n by 1. 

In the program, the estimator is calculated using a vari- 
ation of formula (|23[l that does not contain any integrals. 
Before specifying this variation, let us introduce the ab- 
breviation Njk = (j\N mn \k) and the coefficient Sf m which 
is 1 for I < k < m and otherwise. Then, the following 
expression is equal to Ij23(l for all Hamiltonians Ijl3(l and 
all sets of measurement operators (proof is omitted): 

p n_1 
9 = P , F ' F n = Yl a jj(2k) b Kn-i-k) (27) 

K — U 

with 

hl =(diy T{k+1 2 )T{l + b ^ 

and the recursive relation 

4r 0) = o, 4^ = 

i i 

(n>2) Y^V^/efc at Ar * O-l) 

p=0 q=0 

+8 k 1{ 2 n - 3) i(N j{1 - p) Nf q - N jp N? {1 _ q) )a^ k \ 

+^{2n-A) N j(l-p) N l{l-q) a p n q k ^) ' 

Note, that in order to avoid convergence problems caused 
by low precision floating point data types, the program 
uses data types provided by the GMP library [l(| for the 
calculation of the estimator. 
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